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^ The matrix elements of the zero-range o-force and the finite range 

Gogny-type pairing force are compared. The strengths of the (5-interaction 
for rare-earth nuclei are adjusted. Pairing gaps resulting from different 
pairing interactions are compared to experimental ones. 
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The pairing interaction plays a crucial role in the nuclear mean-field 
model calculations. The most frequently used pairing type is the so-called 
^ ■ monopole pairing interaction ((5r=const model). More realistic residual in- 

^ , teractions are those with state-dependent matrix elements, e.g. 6- and 

Gogny forces. Previously the Gogny-type pairing force were studied within 
- conventional Hartree-Fock-Bogoliubov (HFB) method [T] and in fully self- 

rS I consistent Relativistic Hartree-Fock-Bogoliubov theory (RHB) [^j. Here 

' we investigate different kinds of pairing forces using the Nilsson mean-field 

model and BCS approximation. 

Contrary to the finite-range Gogny force, the 5-interaction is the zero- 
range force. The zero-range nature of the pairing interaction tends to over- 
estimate the coupling to continuum states [ij^. This defect doesn't occure 
when using the Gogny force which is, however, cumbersome to handle; it 
involves sophisticated numerical techniques while the (5-force is relatively 
simple for numerical calculations. The effect of the finite range for the 6- 
force can be easily simulated by introducing an energy cutoff which plays 
the role of an additional parameter similarly like in the case of the monopole 
pairing model except that the 5-interaction is less sensitive to the cutoff. 

These seemingly different pairing interactions should produce in fact 
similar results- the coherence length of nucleonic Cooper pair is of the order 
of size of a nucleus, therefore its structure should be insensitive to the de- 
tails of the interaction in the particle-particle channel The aim of this 



(1) 



2 



paper is to show strong correlations and similarity between matrix elements 
of the pairing part of the Gogny force and of the (5-pairing interaction. The 
pairing strengths of the 6 force for rare-earth nuclei are determined from 
empirical odd-even mass differences. Additionally it is shown that the di- 
agonal matrix elements of the ^-pairing interaction fulfill some dependences 
known from the monopole pairing strength analysis, therefore an alterna- 
tive way of adjusting (^-pairing strengths is given. In order to simplify the 
calculations we assume the Nilsson single particle potential as a generator 
of single particle states and wave functions. 

The finite range part of the Gogny force has the form [31 

Vciru) = + BiPa - HiPr - MiP^P^e-^l^/'^' , (1) 

i=l,2 

where Pa and are the spin and isospin exchange operators respectively. 
In all calculations the DIS parameters [7j have been used and it is assumed 
that only the isovector part of the interaction is active, i.e., Pr = 1 (proton- 
proton and neutron- neutron channel). The zero-range 5 interaction reads 

m 

Vsiru) = Vo{l - &i ■ (T2)/4<5(ri2). (2) 

We limit further considerations to isovector pairing only (T=l). Both in- 
teractions can be then written in a similar form 

V{ri2) = V{\ri-r2\){A + BPa), (3) 

where the radial part is given by (see eg. [Hj) 

-r-^l) = E ^A(ri -r^)— — ^n;(n)yA;.(r2). (4) 
A=0 + M 

It is convenient to calculate the matrix elements of the interaction Q in J 
basis j2j: 

VjMj'M' = {{ah)JM\V{n2)\{a'b')J'M'), (5) 

where a = {laja)- In the case of the same kind of nucleons and for the 
pairing channel one obtains 

Vg = {{ab)m\VG{ri2)\{a'b')m) = 5j^,,d,^,jju,di^,i^, ^ (-1)^(25 + 1) 

S=0,1;A 

{2la + l)(2/fe + l)^2ja' + lV'^jb' + lVxir,r'){A + 3(23 - 1)) 

f L A la' \ j L la' A 1 f /a L ^ \ j ^a' ^a' S 1 , . 
V J \ la' la S j \ I 1 Ja } \ h \ Ja' / '^^^ 
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Vs = ((a6)00|VG(ri2)|(a'6')00) = 5j^i,5,^,jJij,K,i^Ar - r')^ 

A 

The allowed integer values for A in the sums are \la — la' \ 1^ ^ 1^ la + la' ■ 

The simplest pairing model, i.e., the monopole pairing model, assumes 
that all of the matrix elements © in the pairing hamiltonian are equal and 
state independent. 

In order to simplify our studies the calculations were performed in Nils- 
son model (Hj. Formula © consists of two spin terms (5 = and S = \). 
Numerical calculations have shown that the S = 1 term contributes less 
than 0.1% of the total value and can be neglected in further calculations. 

Figure ^ shows diagonal matrix elements {gaa = (aa|^G|5'^'^)) of the 
Gogny force and of the 5 interaction for protons in ^^^Nd, in cases of equi- 
librium and zero deformations. In Figure |2l non-diagonal matrix elements 
are shown for neutrons in ^^'^Nd. We can obtain approximately the same 
values of matrix elements for both interactions considered here by a proper 
choice of the 5 force strength (Vq) which was not done in case of results 
shown in Figures ^ and [21 for the transparency. Strong correlations between 
diagonal and non-diagonal matrix elements of both interactions are common 
for all considered nuclei, either for protons or for neutrons. 

Usually the pairing strength is determined from odd-even mass differ- 
ences for separate nuclei. In Ref. [THj an explicit formula expressing the 
monopole pairing interaction strength g by the average gap parameter is 
derived and a simple dependence of g on the neutron and proton numbers is 
found. In the case of state-dependent pairing interactions the diagonal ma- 
trix elements follow the same particle number dependence as the monopole 
pairing strength. 

Fi gure |3] shows the dependence of diagonal matrix elements gaa 

of the 

(5-force on the particle number 2n equal to the doubled number of levels n 
(the possible number of particles occupying all of the levels up to the level 
n) for the spectrum ^^^Ce. The dependence is of the type 

gaaC^n) = const{hujQg) / {2nf^^ , (8) 

where HuJoe = 41.0[1 ± {N - Z)/A]/A^/^MeV. The plus sign holds here 
for neutrons (/ = N) and the minus for protons (/ = Z). The constant 
values found after the numerical fit of the (5-pairing strength to the mass 
differences are 0.310 Uloon for neutrons and 0.300 faooz for protons. This 
behavior resembles the results of Ref. |10j where it was shown that the 
monopole constant gj\f = 0.284fiujQj\f and gz = 0.290^woz- This leads to 5% 
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Fig. 1. Diagonal matrix elements gaa = i<i<i\VG\5\o.a) v.s. single particle energies 
for protons in ^^^Nd. Vertical lines represent the single particle spectrum. The (a) 
part of the figure corresponds to the equilibrium deformation while the part (b) to 
the spherical case. 



differences in the matrix elements. This may serve as an alternative way of 
adjusting an approximatre (^-pairing strength and may be especially useful 
in case of nuclei far from the beta stability line, where no experimental 
masses are known. 

Pairing gaps measure the lowest excitation energies of even-even nuclear 
systems and determine the binding energy of the coupled pair of nucleons. 
There are a few different ways of evaluating the pairing indicators from 
experimental nuclear masses. The three-point A^^^ indicator was proved 
to be a proper estimate of the pairing gap parameter jTJ. It is defined as 
follows 

A(3)(A/-) = ^[B{M-1) + B{M+1)- 2B{Ar)], (9) 
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Fig. 2. Same as in Fig. ^but for non-diagonal matrix elements which couple to the 
Fermi level F in the case of neutrons in ^^°Nd 



where tt^v" = (—1) is the parity number and B{M) is the binding energy of 
a system consisting of TV particles. Theoretical pairing gaps are calculated 
in BCS model. State dependent BCS equations for the pairing gaps and the 
occupation probabilities are 

The set of BCS equations is solved numerically using a suitable iteration 
procedure. The strengths Vq of the 5 force for protons and neutrons were 
adjusted for the pairing window containing 2^/T5Z{N) levels respectively. 
The values of Vq were determined from the condition of the equality of the 
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Fig. 3. Diagonal matrix elements gaa of the S interaction as a function of the level 
number a for protons (a) and neutrons (b) in ^**Ce. The solid lines represent the 
fit of the matrix elements according to Eq. |S1 



experimental pairing indicator A^^^ and the lowest quasiparticle energy for 
each nucleus and then the values of Vq were averaged for the whole region. 
The strengths of the J-pairing force for rare-earth nuclei obtained this way 
are 0: y^™*""'' = 240MeV • fm^, yj^e"^™"^ = 230MeV • fm^. 

Figure |3 shows the pairing gaps for the 5 interaction (opened circles) 
and the monopole pairing interaction (filled circles) in comparison to ex- 
perimental pairing gaps (Eq|^ for Cerium isotopes. The monopole pairing 
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Neutron pairing gaps calculated with the (5-pairing interaction (opened 
and in the monopole pairing model (filled circles) in comparison to 
indicators (squares). 
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strengths used in the calculations are those of Ref. |^. The state de- 
pendent pairing results in theorethical pairing gaps which differ from A^^^ 
indicators on 10 % while the monopole pairing gaps differ on 20 %. 

In conclusion, both the zero range 5 force and the finite range part of 
the Gogny force have been investigated. It was shown that after proper 
renormalization both forces give the same matrix elements therefore they 
produce the same results. Thus the use of the 6 force is justified and will 
reduce numerical efforts in various studies. 

The 2n~^/^ dependence of the diagonal matrix elements of the (5-pairing 
force can be used in adjusting an approximnate ^-pairing strengths, s 
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